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(57) Abstract 

The present invention 
reduces noise in digital photo- 
graphic images based on the 
assumption that images may be 
decomposed into two types of 
regions, smooth regions and 
edge regions. Smooth regions 
are areas of the image lacking 
any sharp detail, such as blue 
sky. Edge regions are regions 
containing sharp detail, such 
as edges and textured regions 
(such as grass). The present 
method reduces noise in the 
smooth regions by a mathemat- 
ical blurring technique based 
on least squares regression. 
The blurring does not degrade 
the sharpness of the image, be- 
cause there are no sharp details 
in the smooth regions. Edge re- 
gions are left undisturbed to 
maintain sharpness, but the 
noise is less noticeable in those 
regions than in the smooth re- 
gions. The method operates upon the luminance and chrominance component of a digital image. 
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_1_ 

DIGITAL IMAGE NOISE REDUCTION OF LUMINANCE AND 
CHROMINANCE BASED ON OVERLAPPING 
PLANAR APPROXIMATION 

5 

Cross-Ref erence To Related Application 

The present application is related to U.S. 
Patent Application Serial No. fileH on even 

10 date with the present application and entitled 

"A Computer Based Digital Image Noise Reduction Method 
Based on Overlapping Planar Approximation" by 
. M. Kaplan et al . The related application bears the 
Assignee's Docket No. 54,891. 

15 

Technical Field of the Invention 

The present invention relates to an image 
data processing method and, more particularly, to a 
20 method for reducing noise in digital photographic 
images based on smoothing, by a blurring technique, 
the smooth regions of a digital image while 
maintaining the edge regions undisturbed. 

25 BACKGROUND OF THE INVENTION 

To transform a digital image to a finished 
print of high quality requires that noise components, 
introduced because of the transformation of the image 

30 from an image medium to electrical signals be reduced 
and/or eliminated such that the noise does not become 
discernible to the human eye. 

One of the methods receiving widespread use 
in the prior art is related to smoothing the 

35 differences between the values of the gray levels of 
pixels located in neighborhoods. A difficulty 
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associated with this smoothing process is that it not 
only removes the noise components, but it also blurs 
the edges. The edges exist when there is a transition 
region such as a region containing sharp detail, for 
5 example, grass which may define or outline a 
neighboring region which may be either smooth or 

sharply detailed. 

A patent of interest for its teaching of 
noise reduction in a data processed image is U.S. 
10 Patent No. 4.734.770 entitled "Image Data Processing 
Method and Device Therefor" by I. Matsuba. The method 
of that patent treats the first image data set, 
consisting of a plurality of pixels having values 
corresponding to their gray level some of which 
.15 contain noise components, by transforming them into 
second image data to reduce the noise components 
appearing in the second image data. The method 
operates by selecting pixels to be treated and 
obtaining proposed values (candidate values) which are 
20 given to the object pixels based on a stated 

relationship. A probability function is used to 
determine whether a pixel in question should be raised 
to the candidate value or should be maintained at its 
present value.. The probability formula includes the 
25 image energy component as one of its terms .A 

recognition is made for edges of patterns based on 
energy level comparisons, wherein a high energy level 
difference between a central and an adjoining pixel 
would have a high probability of defining an edge. 
30 In the present invention a solution to the 

edge blurring, when performing the smoothing 
operation, approaches the. problem by determining edge 
regions and by leaving those regions undisturbed or 
unoperated upon and by applying a smoothing function 
35 to the smooth regions so as to reduce noise in those 
areas. The method has two responsibilities, one to 
identify whether a pixel belongs to an edge or to a 
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smooth region and. two, to compute a smoothing code 
value for the pixel if it belongs to a smooth region. 
Within this application, these two operations are 
called "segmentation*' and "smoothing," respectively. 
5 "Both of the aforementioned tasks are accomplished by a 
least squares fit of a plane to the image, in a 
neighborhood around the pixel of interest- If the 
goodness of the fit is small, then in that 
neighborhood the image is well approximated by a 
10 plane, and thus must be a smooth region. Further, the 
fit provides a smooth estimate for the code value at 
the pixel of interest which is the value of the fitted 
plane at the pixel of interest. 

15 SUMMARY OF THE INVENTION 

The present invention method reduces noise 
Cgraininess) in digital photographic images by 
operating on the luminance component of an image and 

20 the chrominance components. Generally speaking, this 
is accomplished by a mathematical generalization from 
a scalar field to a vector field, with special 
consideration being given to the characteristics of 
the human visual system. 

25 One preferred method of the present invention 

for reducing digital image noise is comprised of the 
steps of forming a signal matrix of pixel values for 
each primary color of an original image; determining 
if a pixel belongs to an edge region and if not, 

30 computing a smoothing code value from the pixel so as 
to reduce image noise. 

To determine if a pixel of interest belongs 
to an edge region a neighborhood is selected around 
the pixel and a least squares fit is performed to 

35 provide a goodness of fit factor. Next, another least 
squares fitting plane operation is performed using a 
matrix to provide an estimated value of the pixel of 
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interest. The above is continued for all 
neighborhoods around the pixel. All of the estimated 
values are combined with weights based on the goodness 
of fit to provide a pixel value. The operation is 

5 continued for all pixels in the original image and the 
pixel values derived therefrom are used to form an 
output image with reduced digital noise. 

The image is a mathematical field because it 
is a function of spatial variables, namely the x and y 
10 coordinates of the pixel position in the image. The 
luminance image is a scalar field, because it is a 
function having a single value (luminance) at each 
spatial position (pixel). The full 3-color image 
(luminance and chrominance) is a vector field, because 

15 it has 3 values (red, green, blue) at each spatial 
position. To provide a mathematical generalization 
from a scalar field to a vector field, the scalar 
equations need to be rewritten in the proper form for 
a vector field. 

20 From the foregoing, it can be seen that it is 

a primary object of the present invention to provide 
an improved method for removing noise components from 
the luminance and chrominance of digital image signals 
without degrading edge values. 

25 It is a further object of the present 

invention to provide an improved smoothing of the 
pixel values of a digital image in smooth regions 
while not disturbing edge regions. 

It is an additional object of the present 

30 invention to provide a method for reducing noise in 
the smooth regions of a digital image by a blurring 
technique based on a least squares regression applied 
to luminance and chrominance signals. 

These and other objects of the present 

35 invention will become more apparent from the following 
description when takien in conjunction with the 
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attached drawings wherein like characters indicate 
like parts and which drawings form a part of the 
present description . 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a line diagram illustrating the 
smoothing of a near edge pixel by a convolution 
technique. 

10 Figure 2 is a line diagram illustrating the 

smoothing of a near edge pixel by a least squares 
technique . 

Figure 3 illustrates a square neighborhood of 
pixel elements with the target pixel being near an 
15 edge. 

Figure 4 illustrates a matrix of pixels with 
the smoothed pixels being distanced from an edge on 
both sides so as to create a channel of unsmoothed 
pixels surrounding each edge. 
20 Figure 5 illustrates four different blocks of 

(5 X 5) neighboring pixels about one target pixel 
(four quadrant neighborhoods). 

Figure 6 illustrates a matrix of pixel 
elements representing selected neighborhoods for two 
25 target pixels on opposite sides of an edge. 

Figure 7 represents an overlapping of 
corresponding neighborhoods for two adjacent target 
pixels. 

Figures 8A and 8B illustrate a flow chart 
30 depicting the method steps of the present invention. 

Figure 9 is a verbal flow chart verbally 
depicting the method steps of the flow chart of 
Figures 8A and SB. 



35 
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BTsf^CRIPTinN OF THE PREFERRED EM1. :DIMENI 

As previously mentioned, the method of the 
present invention has two major chores. One is to 

5 identify whether a target pixel belongs to an edge or 
a smooth region and the other is to compute a smoothed 
code value for the pixel if it belongs to the smooth 
region. The smoothing process can be done by 
convolution with a smoothing kernel, as is common 

10 practice in digital signal processing. But doing the 
smoothing by least B(iuarcs provides an advantage. 
Referring to Figure 1, wherein the smoothing of a 
pixel by a convolution is shown requiring the kernel 
"neighborhood" to be symmetric and centered on the 

15 target pixel, to provide an unbiased estimate for the 
smoothed target pixel. When the target pixel 
approaches an edge, the neighborhood surrounding the 
target pixel will project pixel miembers onto and past 
the edge surface. A least squares method has no such 

20 requirement. 

• Let Nj^^ a ''^^ * neighborhood of type a, 



centered at pixel location (k,l). 



25 



• Let the 3-component vector Z^j be the color values 

(red, green, blue vector) and Z^- be the 3x3 

noise covariance matrix at pixel (i.j)- 

in • Let X?, be the goodness of fit for the linear 

least squares regression of a hyper-plane on the 

■♦kl.a 

neighborhood Nj^.^^^^ and let Z^j be the estimate 

of the regression at pixel location (i,j), where 

35 pixel (i.j) is some pixel in the neighborhood 

N, , .A multi-variate linear least squares 
kl ,a 
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regression (which is a standard mathematical 

technique) of a hyper-plane finds the vector 

coefficients a, c which minimize the fitting 
2 

error X^, ^ defined by: 
5 JS.J. , a 

(1) 

and the estimated code value at pixel (i.j) is: 

= ^kl + bu i + cuj ( 2 ) 

Compute the noise reduced code value at pixel 
(i.j) as : 



10 



15 



20 



30 



_ , kl.o 



fcU ^ ^ (3) 
where the sums are all kl.a such that pixel 
(hj) € iVi/.Q . 

• Apply the weighting function, 



•■kl o . 



25 "° (4) 

where is an empirically chosen weight for 
each neighborhood type, ii^^ is the number of 
pixels in the neighborhoods is the degrees 
of freedom for the least squares fit (which is the 
number of pixels in the neighborhood a minues 9, 

because the vector coefficients a, S, c have a 
total of 9 components) and is an empirically 
chosen rolloff parameter. 

35 
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To further generalize the methcd, equations 

1-4 are rewritten. Two least squares fits must be 

performed, once for segmentation and again for 

2 

smoothing, so we must minimize two X 's. 

5 

= ^^E^^ ^ K - + i + cf, ;)] ' < - (a-f, + if/ i + cf,;)] 

(5) 

(6) 



(i,y)€A«,. 



Where the matrices M?., M?. are related to the noise 

15 

covariance matrix as explained later. (Note: in the 
notation, E stands for "edge," and H stands for "human 
visual system, ". representing the segmentation and 
smoothing step'^ , respectively.) The weights are 
2Q computed based on the segmentation fit, 

"a 

(7) 

25 

and the least squares estimated code values are based 
on the smc?C>thing fit, 

30 «8) 

The final result is the noise-reduced image code 
value, 



35 
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The most likely application of this further 
generalization would be to define the two weighting 
matrices as follows: 



10 



K^j^iu'EDEUE^K' {i^eDeUe) 



(10) 



(11) 



where Uj,,Ug are unitary matrices. I^E^I^g 

diagonal matrices, and Z^j is the noise covariance 

15 matrix. The unitary matrices Ug,Ug perform a 

rotation in vector space, and the diagonal matrices 
Dg,Dg provide different weights along the 
different axes in the rotated vector space. Thus, a 
likely choice would be to select a Ug which rotates 

20 into a luminance type space (placing one vector space 
axis along luminance Y, and another vector space axis 
along the more important of the two chrominance 
directions), and Dg which weights the luminance most 
heavily, and weights the less important of the two 

25 chrominance components least, ^e'^e selected 
to rotate into a space which is good for locating 
edges, with appropriate weights, which may be a 
luminance space, a principal components space, or some 
other space. Note, that in a principal component 

30 space, by definition U diagonalizes the covariance 
matrix I^j, which simplifies the computation. 
Also, if Uj,,Dj. are chosen identical to Ug.Dgi 
then equations 5 and 6 become identical and the 
computational burden is reduced. If I5j,,Dg are 

35 both the identity matrix, then equations 5—10 reduce 
back to equations 1-4, regardless of the values of 
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Ug.Ug. Thus, rotating the vector space 
accomplishes nothing, "unless the different axes are 
weighted differently. 

Figure 2 illustrates the advantage of 
5 utilizing a least squares technique. The target pixel 
can be located with an off-center neighborhood, 
thereby permitting the target pixel to approach closer 
to an edge region. As Figure 1 shows, a convolution 
' smoothing tends to blur edges, because they "leak" 
10 into neighboring pixels. A least squares method can 
avoid blurring edges by choosing neighborhoods which 
• avoid the edges . 

Referring to Figure 3. an algorithm based on 
the above takes a 5 x 5 neighborhood centered on each 
15 target pixel and performs a least squares fit and. if 
the value (X^ is a sum of squared errors 
normalized for the standard deviation at each pixel) 
is below a threshold value, the algorithm replaces the 
target pixel's code value by the least squares 
20 estimate. This approach is acceptable in the interior 
of a smooth region but fails for target pixels which 
are near the border of a smooth region. This failure 
is caused by the 5x5 neighborhood overlapping the 
edge or border of the smooth region, which in turn 
25 results in a large X^ that exceeds the threshold 
value and therefore the pixel of interest is not 
replaced with the pixel code value or in other words 
is not smoothed. The result of continuing that 
process around an edge is shown in Figure 4 such that 
30 a channel or roadway of unsmoothed pixels around an 
edge line will result. 

An improvement is achieved by considering 
several neighborhoods for each target pixel and 
selecting one of the neighborhoods which results in 
35 the most preferred fit. For example. Figures 5A 
through 5D illustrate a target pixel being a corner 
pixel for four 5x5 neighborhoods. These 
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neighborhoods can be called quadrant neighborhoods and 
are used to determine the value of which has the 
smallest value. This value can then be compared 
against the threshold value to determine if the pixel 

5 value should be replaced with a code value. If none 
of the neighborhood values is below the threshold then 
it is assumed that the pixel is truly an edge pixel 
and a smoothing is not applied to it. If at least one 
neighborhood value is below the threshold, then the 

10 neighborhood with the lowest value is used to compute 

the smoothed value. 

In Figure 6, the selected target pixel 
neighborhoods, for two target pixels, are shown 
positioned on opposite sides of an edge. The result 
15 of this is that for pixels near an edge, the pixels 
will be smoothed by incorporating them into 
neighborhoods entirely on the correct side of the 
edge. This method provides proper noise reduction 
right up to the edges (as opposed to leaving bands, of 
20 noise surrounding the edges), but it also produces an 
artifact which makes the picture looK "blocky". The 
artifact is caused because the algorithm has a 
tendency to use one type of quadrant neighborhood 
"such as the upper left" for several adjacent target 
25 pixels, and then to switch to another type of quadrant 
neighborhood (such as a lower right) for. the next 
several target pixels, and so on. 

To eliminate the "blocky" look, all of the 
neighborhoods containing the target pixel were used. 
30 For 5 X 5 neighborhoods, for each target pixel there 
are 25 neighborhoods containing the target pixel. 
This is intended to provide a more continuous 
behavior, to reduce blocking, because of the larger 
number of neighborhoods used. 
35 unfortunately, the "blocking" persists even 

after using all neighborhoods . This was attributed to 
the fact that for two adjacent target pixels, there 
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are two sets of 25 neighborhoods (one set for each 
target pixel), and each neighborhood in one set 
contains mostly the same pixels as several 
neighborhoods in the other set (for example see 
5 Figure 7). Thus the two sets of least squares fits 
produce similar results, and when the algorithm 
selects the neighborhood with the lowest X for each 
target pixel, it tends to pick nearly identical 
neighborhoods from each of the two sets of 

10 neighborhoods. As the algorithm scans the picture 
selecting new target pixels, it continues to select 
similar neighborhoods. Occasionally the selection of 
neighborhoods flips to a radically different choice, 
and then this new selection persists over may target 

15 pixels. A noticeable blocking artifact is produced 

where the selection flips. This is basically the same 
blocking artifact as described above for quadrant 
neighborhoods. 

Since the blocking artifact appears to be 

20 caused by the instability in the selection of the 

single best neighborhood for each target pixel, the 

algorithm was modified to avoid the need to make a 

selection. Instead of selecting one neighborhood, the 

algorithm uses all neighborhoods of a target pixel, 

25 weighting them according to "goodness," applying large 

2 

weight for small X (good fits) and small weight for 

2 

large X (poor fits). This eliminates the blocking 



artifact 



2 

The weighting function w(a,X ) (equation 



2 

30 ^) originally based on the X probability 

distribution, which is statistically appropriate when 

the model surface (a plane) correctly describes the 

2 

data surface (image), with the fitting error X 
resulting only from noise and errors in choosing the 

35 fitting parameters. However, analysis showed that the 

2 

fitting errors X were statistically distributed 

2 

very differently from the X distribution, even in 



10 



15 
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the smooth regions. This indicates that the model 
sijrface (plane) is not an adequate description of the 
data surface (image). This is not surprising, since 
the data surface has a very curved and irregular 
5 shape, even within the smooth regions. To deal with 
■ this difficulty, the M weighting function of 
equation 4 was substituted for the original X 
distribution based weighting function. 

neighborhood size selection is a remaining 
factor. The larger the neighborhood, the larger the 
noise reduction, because many pixels contribute to the 
smoothing. However, a large neighborhood size is 
unable to deal with details in the picture smaller 
than the neighborhood size. This is because for small 
details there is no neighborhood which fits entirely 
within the detail, and thus the algorithm considers 
the detail to be an edge region, and provides no 
smoothing. This difficulty is dealt with by using 
neighborhoods of several sizes. The algorithm weights 
different sized neighborhoods differently, with the 
result being that large smooth regions of the picture 
are smoothed primarily by large neighborhoods. 
Smaller smooth regions are smoothed primarily by 
smaller neighborhoods, because the large neighborhoods 
don't fit within the region, thus have large X and 

small weight. 

Finally, it should be noted that the 
algorithm includes neighborhoods of size 1x1. 
Although this can be mathematically described simply 
as 1 X 1 neighborhoods, the intent and result are 
somewhat different than for larger neighborhoods — 
the 1 X 1 neighborhoods provide the mechanism by which 
edge regions are permitted to remain unsmoothed. In 
an edge region, all neighborhoods larger than 1 x 1 
yield large when the least squares fit is 
performed. Strictly speaking, it is not possible to 
fit a 1 X 1 neighborhood with a plane, because there 



20 



25 



30 



35 
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is only one 3-component data point and nine fitting 

parameters (making the fit mathematically 

underdetermined) . Instead, a perfect fit is assumed 

(i.e. the fit estimate of the code value at the single 

5 pixel is exactly equal to the original code value, and 

= 0). The result is that in edge regions, only 

2 

the 1x1 neighborhoods have small X and are 
heavily weighted, and they pass the original picture 
information unsmoothed. 
10 Referring now to Figures 8A, 8B and 9 in 

conjunction, from a multi-band original digital image 
(Blue, Green and Red) a target pixel is selected. For 
each target pixel, a neighborhood is selected. The 
total number of neighborhoods for the preferred 

15 embodiment is 588. For each neighborhood selection 
the least squares best fitting plane is computed using 
the "Edge" matrix (equation 5) to yield the quality of 
fit X^. Additionally, for each neighborhood 
selection the least squares best fitting plan is 

20 computed using the "Human Visual System" matrix 

(equation 6) to yield the estimated value (z) at the 
target pixel. The process is iterated for all 588 
selected neighborhoods. All 588 estimates (z) are 
combined with weights based on the quality of fit 

25 (X^), according to equation 3 to yield the output 
image- pixel value at the target position. This 
process is iterated until all pixels in the original 
image have been processed. 

This approach generalizes the approach of 

30 Assignee's Docket No. 54,891 as follows. 

First, a full noise covariance matrix 
(Z..) is used, and thus it takes into account 
correlations of the noise in the different channels. 
Second, the weights w(a,X^) are computed for all 3 

35 channels together, not independently, and the same 
weight is used in each of the 3 channels. Computing 
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the weights by combining the channels enhances the 
ability to recognize edges and textures. Third, it 
takes into consideration the human visual system. 

It is desirable to further generalize the 
5 equations, to take into consideration the 

characteristics of the human visual system. The 
present method is a two step method. First, it 
segments the image by computing the weights (equations 
5,7). Second, it smooths the image code values 
10 (equations 6,8,9). There is a significant difference 
between these two steps that must be considered. The 
smoothing step computes new image code values, which 
will be directly observed by the human eye, so special 
attention is paid to the luminance channel because of 
15 its greater significance to the human visual system. 

On the other hand, the segmentation step 
represents an effort to locate edges and textures in 
the image, and does not compute new image code values, 
and thus luminance is less relevant in the 
20 segmentation step. If luminance is significant in the 
segmentation step, it is for a more subtle reason. 
The human visual system has evolved into a highly 
efficient system for understanding the visual world, 
and thus perhaps its emphasis on the luminance channel 
25 is a commentary that when trying to understand a 
scene luminance is given extra consideration. 

While there has been shovm what is considered 
to be the preferred embodiment of the invention, it 
will be manifest that many changes and modifications 
30 may be made therein without departing from the 

essential spirit of the invention. It is intended, 
therefore, in the annexed claims, to cover all such 
changes and modifications that may fall within the 
true scope of the invention. 

35 



15 
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CLAIMS: 

1. An overlapping planar approximation 
method for reducing digital image noise in multi-band 
digital images comprising the steps of: 
5 (a) selecting a target pixel (which is 

a vector with a number of components equal to the 
number of bands in the image) from a set of pixels 
representing an image of interest; 

(b) selecting a neighborhood of the 

10 target pixel; 

(c) computing a least squares best fit 
plane of pixel' vectors for the selected neighborhood 
of step b to provide a goodness of fit (X ) and an 
estimated vector value of the target pixel; 

(d) iterating steps b and c for various 
types (sizes) and positions of neighborhoods of the 
target pixel; 

(e) summing all estimated values of the 
target pixel with weights based on the goodness of fit 
(X^) to provide an output image pixel for the target 
position; and 

(f) iterating steps a through e for all 
pixels from the image of interest, so as to form a set 
of reduced noise image pixels. 

25 2. The overlapping planar approximation 

method, according to Claim 1. wherein the goodness of 
fit of step c is defined as: 

"^^ (1) 

and the estimated code value vector of a target pixel 
positioned at (i.j) is: 



20 



35 



z 



•J (2) 
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and the computed noise reduced code value vector of 
the target pixel at position (i,j) is- 

£.^(°.vi,..)^-.;'' 

where the sums are over all kl.a such that pixel 
€ iVi/^a sind where the weighting function is: 

10 

(4) 

15 wherein Nti^^ is a neighborhood of type q centered at 
pixel (^,0, Ho is an empirically chosen weight for 
each neighborhood type, rto is the number of pixels in 
the neighborhood, i/^ is the degrees of freedom for the 
lease squares fit, is an empirically chosen rolloff 

20 parameter of order unity, ^0* is a code value vector 
(having one component per band of the multi-band 
image). Eij is the noise covariance matrix at 
pixel (^i)' xli^a is the goodness of fit for the 
linear least squares regression of a plane on 

25 neighborhood A'^/^q , and r,^''** is the estimate of the 
regression at pixel location , where pixel (^i) 

is some pixel in the neighborhood A'ti.o. 

3. The overlapping planar approximation 
method, according to Claim 2, for the case of the 

30 noise covariance matrix E,j being diagonal to thereby 
partially decouple the equations for the various bands. 
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4. The overlapping planar approximation 
method, according to Claim 2, wherein two least 
squares fits are substituted in place of the single 
least squares fit by replacing equations 1 through 4 
5 by the following equations: 

Xhia = J2 ' + ^ki i + ^ij)] ^ M.f [z,j - (af, + bf, i + cf,i)] 

(5) 

10 

xlkLa = E [^ij - (SS + i + 4 i ) ] ^ - + fcf, i + cS j) ] 

{»,y)€A'»i.ci 

(6) 



15 



20 



and the estimated code value of a target pixel 
positioned at (i,j) is: 



and the computed noise reduced code value of the 
target pixel at position (i,j) is: 



(7) 



where the sums are over all kl.a such that pixel 
30 (».;) e A'iu where the weighting function is: 



35 
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where and M-^ are arbitrary matrices, chosen to 
optimize segmentation and smoothing, respectively. 
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5. The overlapping planar approximation 
method, according to Claim 4, wherein the matrices Af^^ 
and MJI ate defined as follows: 

5 * = {u^EDsUEyKj' (uIDeUe) 

(10) ' 
(11) 

10 • 

where Ue^ are unitary matrices . De.Dh are 
diagonal matrices, and is the noise covariance 
matrix such that the unitary matrices provide rotation 
> of the vector space axes and the diagonal matrices 
15 attach different weight to the different axes in the 
rotated vector spaces. 

6. The overlapping planar approximation, 
method, according to Claim 5, with the unitary 
matrices chosen to rotate into a luminance/chrominance 

20 vector space, and the diagonal matrices chosen to 
assign greater weight to the luminance axis in the 
rotated spaces. 

7. The overlapping planar approximation 
method, , according to Claim 1, wherein said goodness of 

25 fit (X^) is used as a measure of edge detection, 
texture detection, and image segmentation. 

8. The overlapping planar approximation 
method, according to Claim 4, wherein said goodness of 
fit (X^) is used as a measure of edge detection, 

30 texture detection, and image segmentation. 

9. The overlapping planar approximation 
method, according to Claim 5, wherein said goodness of 
fit (X^) is used as a measure of edge detection, 
texture detection, and image segmentation. 

35 
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10 A method of image smoothing and noise 
reduction by least squares fits, comprising the steps 

(a) selecting a target pixel from a set 
5 of pixels representing a multi-band image of interest; 

(b) selecting a neighborhood containing 

the target pixel; 

(c) computing the least squares best 

fit plane for the selected neighborhood of step b; and 

(d) using the least squares fit 
estimated value of step c as the smoothed output value 

for the target pixel. j ^^^ce 

11 The method of image smoothing and noise 

.eduction, according to Claim 10. -^"^^VJ^'/^^f 
sauares fits are substituted in place of the least 
scares best fit of step c by the use of the following 

equations; 

(•j)€Aii.« 

20 (5) 
XW= Z [^-.-(°t, + tS,+ct,i)]\f.?[r,-(at, + 6f,i + ct,^ 



)€.Viti.o 

25 



(6) 



30 



and the estimated code value of a target pixel 
positioned at (i,j) is: 



and the computed noise reduced code value of the 
target pixel at position (i,j) is: 



Ziofa.xiiu) (8) 
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where the sums are over all kl.a such that pixel 
(ij) £ Nti,a and where the weighting function is: 

(9) 

where M,j and are arbitrary matrices, chosen to 
optimize segmentation and smoothing, respectively. 
10 12. The method of image smoothing and noise 

reduction, according to Claim 11, wherein the 
matrices M^j and M,^ are defined as follows: 

15 (10) 

(11) 

20 where I'e^I'h are unitary matrices , D£, Djj are 
diagonal matrices, and E,j is the noise covariance 
matrix such that the unitary matrices provide rotation 
of the vector space axes and the diagonal matrices 
attach different weight to the different axes in the 

25 rotated vector spaces. 

13. The overlapping planar approximation 
method of Claim 1 applied to multi-channel digital 
signals other than multi— band digital images. 

14. The overlapping planar approximation 
30 method, according to Claim 2, wherein the weighting 

function is selected to provide the desired result. 

15. The overlapping planar approximation 
method, according to Claim 2, wherein the selected 
neighborhood sizes of the target pixel are 21 x 21, 

35 
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11 X 11, 5x5, 1x1, the empirical weights are: 
TV21X21 = 64. W„xn = 16. W^sxs = 16. Wi^i = 1, 
and the rolloff parauneters are: 

r21x21 = nixll = »"Sx5 = nxi — 1- 

5 
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